Popular Computing 


The world’s only magazine devoted to the art of computing. 


Volume 9 Number 4 
April 1981 


A Numberprint For Pi 


Pow fle 


[a co FF WW OO Ww On 


Oo Vr Ye Ph 2 ea SG Ss © 


=) 


w 
NT AX oOo. F&F NAN wh Oo WwW A WwW OF MO 06 80 LF Dees omvunn. 


(OR (Yan eo) = ES (oo) No) Noy ONoji =e No SS ey TS! op, RD) 

Wim Ml Oy Corre © Ur NM “GY CO Oro) Snr Oy 1 TON CoC NmCcomiet 
“P56 25M OF. 0 fF NW FY © WS }W BOO Ul Micro 
“3a> FF NN fF Oo COC OC WU NAN @© eA OO OOS Seco wie 
Se tt PrP Cf NY OF ow Sa oe An) 6D CY 9 2 eee ees 
ee rent CN es) te ON to OR, ad 

Nee ©) (OM SN Cone Go Cx CON Onn NOME mn) 

oS NM ONO. F&F WM OY @ fF CY OY LE Go) (Co ee iomeonmncs 


fo) 


An example of a scheme for converting 
a fraction (in this case, 9 digits 

of the mantissa of pi) from decimal 
to binary. 


In each of the Numberprints shown, 
the origin is indicated with an 
arrow. The heavy lines indicate 
legs of the print that are traversed 
more than once. 
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Numberprints 


Suppose we take a number, like pi, and convert it 
from its familiar decimal form to binary. The rules 
for converting integers and fractions are different. In 
Figure S, a scheme for the conversion is shown. The 
fraction is doubled, and whenever the doubling spills over 
across the decimal point, a "1" is generated in the binary 
number, starting at the binary point; otherwise, a "Oo" 
is generated. In the illustration, the first 20 bits of 
pi are developed, which is about all that can be expected 
from a 9-digit approximation to pi. 


We thus arrive at: 
pL = 11.00100100001111110110101010001000... 


Now, 1f we interpret these bits two at a time 
according to this scheme: 


00 North 

Ol East 

10 West 

at at South, 
then we have a means of picturing numbers. The cover of 
this issue shows the Numberprint for pl. The Numberprints 
for e (the natural logarithm base) and the square root of 2 
are also given. 

Flowchart R indicates a possible algorithm for 

generating Numberprints, from an input of 29 decimal places 


of a fraction. From this input, it should be legitimate 
to generate about 50 legs of the Numberprint. 
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Numberprints for e and the square root of 2. 


The pattern is readily generated in any of 


the BASICs, and could be produced dynamically 


on a CRT display. 
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1.00000 

-51763 

-26105 

. 13080 

-06543 

192 03272 

384 .01636 

768 .00818 
1536 00409 
3072 .00204 
6144 .00102 
12288 .00051 
24576 -00025 
49152 .00012 
98304 .00006 
196608 .00003 
393216 -O0001 
7864 32 - 00000 
1572864 -00000 
3145728 -00000 


00000 00000 
60498 
62476 
93734 


81019 


80902 
23844 
62584 
81656 
34632 70924 
23715 
96636 
56100 
98543 
84314 
30767 
67381 
01535 
10073 


22792 
12080 
06125 
53073 
26538 
13269 
56634 
78317 
39158 
19579 
59789 
79894 
39947 
19973 


12208 
12741 
62875 
75408 
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78542 


The table above is furnished by John W. Wrench, Jr., to 
replace the one on the top of page 11 of our issue No. 12. 
N is the number of sides of an inscribed polygon in a 
circle of unit radius; S is the length of a side; P is the 


perimeter. 


Mr. Wrench says "I believe that the second 


part of the table is also infested with errors." 


00000 
29637 
39432 
41427 
27622 
15090 
29103 
81734 
19148 
10109 
89519 
91038 
79334 
30341 
68338 
27853 
92733 
33953 
69258 
03084 
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Ulam-3 


In issues 13 and 19 there was discussion of a 
problem due to S. Ulam: 


In the sequence 1, 2, 3, 4, 6, 8, 11, 13,. 
each new member can be formed in one and only 
one way by adding two earlier numbers of the 
sequence. 


The first 110 terms of this sequence were given, calculated 
by a program of Associate Editor David Babcock's. 


The accompanying flowcharts show a scheme for extending 
the idea to numbers that can be formed in one and only one 
way by adding three earlier numbers. 


The first few terms of this sequence (call it Ulam-3) 
are given here: 


1 2 3 4 6 7 8 
19 al 22 23 38 58 59 
TA 77 91 92 108 124 125 

127 163 178 198 213 230 2h6 


247 248 249 265 266 282 283 


299 316 336 337 352 369 371 
386 404 406 422 


The 28th through 31st terms (underscored) are 
conjectured to be the only instance (other than the first 
four generating terms) of four consecutive numbers. 


The problem lends itself to coding in almost any 
language (even integer BASIC will carry it quite far). 
The suggested solution, shown on the following two pages, 
slows down considerably as the numbers increase. Hence, 
to extend the sequence another 100 terms would require 
either a machine language approach or, more likely, one of 
those analytic solutions that brilliant problem solvers 
like David Ferguson keep coming up with. 
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Knuth 


One recurring theme in our pages has been the notion 
that there are problems for which the computer is the only 
tool for solution. There are indeed such problems (such 
as certain combinatorial problems) but many of our attempts 
to devise good computing problems were thwarted by having 
someone come up with an ingenious analytic solution, thus 
spoiling the fun. It just happened again: 


In issue 95 there appeared a straightforward solution 
to one of the Penny Flipping problems, using an approach 
that lent itself to machine language on a microcomputer. 
This was clearly a computer problem and for most people it 


stili is. Not, however, for the master; Prof. Don Knuth 
writes as follows: 


When PC95 arrived last week I took another look at penny flipping ... and I got lucky! 


Enclosed is an extension of “Table W”. For each 7 J list the period length p(n) when the penny 
flipping starts to repeat, and also any smalicr number of steps where the coins all come out heads 
up but the last move was not to flip the whole stack. These latter “sporadic solutions” are rare, 
aithough they occur for infinitely many 7m; it isn’t hard to prove that you get a sporadic solution 
of length p(n) — 2k —1 whenever n =1-+-2-++----+ 2k = k(2k-+1) and k > 1. All of the other 
sporadic solutions on this list are of the form 4 p(n) — 1; such solutions occur when the pennies in 
the second half of the period are complemented from their state in the first half of the period. It 
seems that there are two different sporadic solutions only when n = 21; but it will almost surely be 


impossible to prove this conjecture rigorously, or to prove that all sporadic solutions are of these 
two types. 


I plan to give this nice problem to our grad students the next time I teach Stanford’s programming 


seminar. It illustrates the way a good idea can speed up computations even for problems that 
appear to have no simple structure. 


In fact, the use of some elementary ideas from arithmetic and group theory make it possible 
to compute Table W as fast as my hardcopy terminal could print it at 30cps, even when I was 
timesharing the computer with many other people. The casc n = 113 required about .6 of a 


second, by comparison with the method of your program which would have taken more than two 
days. 


The program that produced this table was about 150 lines of ALGOL code. As you know, I 
have nothing against assembly language programming—once I had to write some code that was 
to be executed a trillion times, so that saving 100 nanoseconds in the loop meant saving hours 
of expensive computer time—bul, sometimes it helps more to have a good idea. Yes, Euclid’s 
“algorithm pays off, even when flipping pennies. 


It is clear, I think, that Prof. Knuth must have some 
very clever way of solving this problem, but we may never 
know what it is. He mentioned that his solution for 
N = 4096 took 25 minutes of CPU time. The method that 
appeared in issue 95 would require something on the order 
of 78 centuries. 


Anyway, it is gratifying to find that not only has 


Prof. Knuth forgiven me for my blunder in issue 92, but 
is willing to contribute to the magazine once again. 


But I can't resist the temptation to try, just once 
more, to develop a simple problem that must be solved by 
computer. Try this variation of Throwback (the original 
version was on the cover of issue 55): 

Consider all the integers, starting with 3: 
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The number at the head of the stream is called the leader. 
The leader is to be thrown back in the stream the number of 


Thus, the first few moves are these: 
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When will the number 100 first appear as the leader? 


Can this number be found without a computer? 
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121 

240 119 
130 
1176 
300 
1440 
1428 


> §76 


2128 
320 
882 875 440 


: 968 483 


2622 
432 
1900 
780 
4320 
756 
4466 
26400 
11284 
2688 
20790 
28560 
16170 
16848 
1€872 
83980 
13104 
15960 
30668 
5880 
13330 
230384 
62100 
2484 
1410 
264960 
32634 16316 
12900 
27744 
2600 
2226 
1156680 
75900 75889 
11760 


13680 
493696 
14868 
140400 
14640 
135408 
232848 
61072 25535 
1034280 
54648 
67536 
138720 
8694 4346 
260680 
349888 
47620 
91980 
186480 
7182000 
20520 
5544 
28080 28067 
1161774 
336000 
5320728 
1883376 
68475 
1585836 
285090 
1475760 
456228 
2265120 
1226420 
75600 
301938 
380880 
539028 
78960 
1938000 
2428416 
30264 
19208 9603 
67320 
273000 
8243620 
1105272 
§19120 
18258240 
2425500 2425485 
2281968 
38520 
21168 
6992602 
23760 
1210344 
2173248 
88750200 


105336 
426040 
341040 
2069496 
12198368 
68544 
13543200 
49832640 
717360 
127920 
11561760 
7980000 
89208 
748030 
2376192 
143190 
146900 
149340 
91872 
3527160 
7211880 
256397400 
228480 228463 
13941120 
139104 
36140 
399840 
473760 
161880 
2793648 
290304 
36090 
7901520 
141190560 
2051280 
2199240 
198000 
3234420 
139840 
136256904 
17109400 
39776100 
35865648 
414480 
853200 
11514780 
207486720 
41972700 
100440 
665705040 
3247200 
4643100 
1593600 
13916778 
3153024 
2413320 
437920 


805410 805391 
§8480 29239 
4650240 
171216 
4395300 
658240 
$559216 
42207360 
2536072 
4773600 
230796720 
64853344 
192150 
2104118016 
1361600 
5468400 
4319700 
2320871280 
14277025980 
282156840 
7556724 
322560 
1621200 
939736 
2550600 
939624 
74143314 
25945920 
802368 
891480000 
8842392 
52940160 
97310080 
156143232 
6371400 
4766402392 
5738040 
29213184 
60969480 
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232: 
Beh) 9 
234; 
2356: 
236: 
237: 
238: 
239: 
240: 
241: 
242: 
243: 
244; 
245: 
246: 
247: 
248: 
249: 
250: 
abies 
252: 
eis} 0 
254: 
Zone 
256: 
§12: 
1024: 
2048: 
4096: 


438918480 438918459 


33359100 
21452704 


508725458 304362728 


179439000 
9876240 
6188832 
3098760 
364496 
7015008 
853723200 
704888340 
6717720 
187320 


: 963217920 


376200 
§9018544 
26164020 


299592 149795 
6718860 
21307200 
255633840 
§554080 
147475020 
1432080 
3584220 
4663360 
8859060 
9611868 
44568720 
35006400 
6632170 
401329242600 
763408800 
404796 
2603370000 
6125400 3062699 
3230760 
184675680 
6294720 
3077000 
13599180 
589680 


999827664 999827641 


19431000 
44553600 
738017280 
571936768 
43231518720 
964767498240 
1575568834560 


Knuth's results for the 
Penny Flipping II problem 


Fred Gruenberger 


Swan Song 


This is the last issue of POPULAR COMPUTING in its 
present form. The title will now be used by McGraw-Hill; 
the first issue with the new format by that group will 
emerge in October, We certainly wish them well. 


I would like to summarize the adventures of publish- 
ing a highly specialized magazine for eight years. 


What we had in mind when we started was a small 
business (actually, a self-supporting hobby) that would be 
run right; no customer would ever have legitimate grounds 
for complaint. Since the product is a magazine, the 
business is quite different from other businesses; for 
example, a magazine is the one product that people custom- 
arily pay for long in advance. We guess that, since we 
were dealing mainly with computer people, we could safely 
eliminate nonsense like "Renew now and get this FREE 
booklet." There is, of course, no way to tell whether 
or not this approach was the right one, but we felt that 
we had a fairly exact count of just how many people in the 
world were interested in our view of computing. 


We ran no paid ads, and did little or no advertising 
ourselves. We figured that the word would spread by 
osmosis and indeed it did. A mention in Martin Gardner's 
column in Scientific American was a big help. 


The main idea was to produce the kind of reading 
matter each month that I would like to read myself. Thus, 
"Letters to the editor" saying how wonderful the magazine 
is seemed inappropriate, space-wasting, self-serving, 
and--perhaps-~redundant. "Articles that were thinly 
disguised sales pitches were rejected, as were earth- 
shaking discoveries that computers can perform address 
modification. 


We had our policies and our goals--the latter were 
spelied out in issue number 30--and they are still valid. 
Some of them were: 


To encourage computing for fun. 

To demonstrate what to compute as well as 
how to compute. 

To emphasize the testing of programs. 

To show that there is always a better way. 

To call attention to the pitfalls and booby 
traps of our machines and languages. 

To promote the art of computing. 

To offer new problems for computer solution. 


That last category will provide amusement for many years, 
as our original problems get picked up and reappear in other 
magazines and books. 


We had to pick a level to aim at, somewhat above those 
to whom indexing and looping is a new and complicated notion, 
but below the level of the professionals who are creating 
our new systems and languages; it works out to about one 
semester of calculus. The mail was about equally divided 
between those who claimed they couldn't understand a word of 
it and those who objected to paying for such obvious trivia. 
In any event, we decided to leave it to others to teach 
elementary programming, rudimentary problem solving, and 
all dated information. With some exceptions, all the 
material that has appeared in the magazine has been timeless. 


When producing a magazine more or less single-handed 
(in the first 96 issues, just 26% of the material was 
contributed by others), there are never any nagging worries 
about "what shall we do today?" As each issue is hauled 
to the postoffice, the pressure is already on to get 
eracking with the next issue. With this final issue, 
there have now appeared 1784 pages, of which 1320 pages 
were made up of original material written by me. Nearly 
300 new computing problems were presented, with solutions 
for most of them. The names of some of these problems 
commemorate our travels (e.g., the Repulse Bay Trip, in 
issue 49). 


The bragging above should not detract from the 
amazing contributions of others. David Babcock's 
brilliance influenced nearly every issue. We were proud 
to present the advanced thinking and competent writing of 
Richard Hamming, Thomas Parkin, and Daniel McCracken. 

Bill McGee's book reviews were a joy; he never tried to 
avoid the books that were difficult to review. 


And, of course, the whole venture would have been 
impossible without the support, encouragement, and darn 
hard work of the cute publisher, whose idea it was in 
the first place. 


When the magazine was started in 1973, many people 
took pains to convince us that we would run out of 
material within six months. In the eight years, there 
was exactly one month when we had to scramble for material. 
Most months there was a steady backlog of 40 to 60 pages 
of good stuff ready to go. 


We had two slogans: 


"The way to learn computing is to compute" 


"The world's only magazine devoted to the art 
of computing" 


For the latter slogan, we managed to devise some 70 distinct 
permutations (e.g., "The magazine for people who enjoy 
computing") and proceeded to run through them systematically 
on the covers. No one seemed to notice. 


The magazine developed its ultimate character after 
about 12 issues. There were many appeals to change its 
character into something else; the most intense of these 
efforts was the bombardment by the APL fanatics, although 
in the last two years that has vanished (now it's the 
Pascal fanatics who are busy). I venture to predict that 
by 1982 another magic language will come along and we will 
hear all the old claims repeated again. Am I the only 
one who recalls going through the same charade with PL/I, 
FORTH, etc., etc.? 


It took quite a while to get the printing problems 
licked. Issue number 6, for instance, was delayed when 
we found that the printer who held the masters had calmly 
closed his doors and gone fishing for a week. Issue 
number 66 has given everyone fits, since it was trimmed 
1/4 inch larger in both dimensions than all other issues. 


We tried our best to keep the literacy level up 
(you may have noticed that many articles in the current 
personal computing periodicals seem to be written by 
junior high school dropouts and show no evidence of any 
editing) and insure a decent standard of clear English. 
The keen eyes of our Associate Editors saved us from 
many a glitch. 


After a magazine gets listed in standard reference 
works, its name gets on other lists, such as mailing lists, 
and these breed like rabbits. Not only do press releases 
arrive by the bale, but also regular issues of the Paint 
and Hardware Journal and ads for wire wrap machines, what- 
ever they are. And among this detritus were the weekly 
letters that said "My teacher told us to write a paper on 
computers so please send all the information and hurry."-- 
only those letters weren't spelled that well. 


A magazine subscription list is an amazing thing to 
manage. Think about it: it is, by definition, highly 
volatile and full of unusual twists. Subscriptions that 
arrive from many foreign countries are widely separated in 
both time and space from the check that pays for them. 

Most company subscriptions arrive through an agency, and 
these subscriptions account for about 95% of all circulation 
problems. The communication channels in most large 
corporations are so tangled that the person who wants the 
magazine and requests it, and the agency clerk who relays 
the order, are poles apart. On one infamous subscription, 
we accumulated over an inch of correspondence getting it 
squared away. 


We will not miss the monthly interface with the 
Postal Service (sic). Anyone who deals regularly with 
its inmates can recount tales of horror. When we began, 
in 1973, we paid 4.8¢ per copy and they did all the work; 
when we ended eight years later it was 8.4¢ per copy and we 
did all the work, under a set of rules you wouldn't believe. 
Each month the cost goes up in one way or another, the 
so-called service attenuates even more, and the rules of 
their game become ever more weird. One example will 
suffice. In our early years, we printed our envelopes 
with standard messages about "Return Postage Guaranteed" 
and "Change of Address Form Requested." We gave that 
up. In the last few years, when the address on a copy 
did not precisely match the occupant there, the postoffice 
charged us 28¢ to return the copy. and then only with the 
fascinating message "undelivered." What one wants, and 
what used to be supplied for 10¢, was the new address. 


In our first two years, we received two flattering 
reviews in Computing Reviews, a publication of the 
Association for Computing Machinery. By dint of much 
howling and screaming, we managed to coerce them into a 
third review in 1978. Considering that POPULAR COMPUTING 
ran articles by two national presidents of ACM (neither 
of which they reviewed), not to mention the two marvelous 
articles by Don Knuth, one would think that a subsidized 
journal that purports to review the field would be able 
to manage more than three reviews in eight years. 


The time frame for our magazine overlaps the deluge 
of personal computers and the emergence of some 30 or so 
periodicals that attempt to serve this new market. Along 
the way, several dozen magazines were born, got to issue 
2 or 3, and quietly died. Even among those that have 
survived, none are edited by seasoned computer people. 

The various editors display astonishment (and, frequently, 
horror) at discovering that the number .1 just does not 
exist in binary; that programs can be and are copied; 

that computers can and do break down; that there is no 
upper bound to the number of BASIC renumbering programs 
that the world will buy; that unscrupulous people will 
try to pass off other people's work as their own; that 


2 
(sine x) + (cosine x) 


adds to unity for lots of values of x; that interactive 
programs to guess a number between 6 and 8 can be sold 
as teaching tools" over and over. It seems that 

snake oil outsells medicine 100 to 1, as it always has. 


To paraphrase Norman Sanders, it would be churlish 
of me not to acknowledge the contributions of California 
State University, Northridge, which provided me with many 
long committee meetings at which I could appear to be 
alert and attentive the while devising new computer 
problems. 


A footnote to history: we determined, after 
exhaustive research, that Eagle "Prismacolor" Light 
Blue #904 pencils will disappear from Itek cameras 
and most photocopying machines better than any 
other such pencil. 


Since there are always many folders lying around 
labelled "Work in progress," there is still a lot of 
unfinished business, as for example: 


1. When and if Contributing Editor Thomas R. Parkin 
generates Part III of his definitive essay on Time (the 
first two parts were in issues 63 and 86), it may go 
unseen and enheralded, except by Tom and me. The main 
item left dangling is a magic formula that will convert 
dates from various irregular calendars (such as our present 
calendar) to the orderly scheme of one Joseph Scaliger. 
Many people, including competent astronomers, have tried 
to devise such a formula and failed, but Parkin cheerfully 
maintains that he can do it. 


@. There was going to be an article about the 
unique features of the IBM 1620, and particularly its 
ingenious method of doing arithmetic by table lookup. 


3. The work on placing isosceles triangles on a 
lattice (in issue 96) was intended to lead to a nice 
difficult computing problem. 


4, Then there is the Foreign Exchange Problem. The 
daily paper lists the value of some 30 foreign currencies 
in terms of U.S. dollars (as on the cover of our issue 
number 83). For example, the rates for Peru and Israel 
have been: 


Peru Israel 


May 8, 1980 270.24 

The number 
Aug 13, 1980 291.30 

of units of 
Sep 10, 1980 293.94 

the other 
Nov 7, 1980 317.42 

currency that 
Jan 8, 1981 334.35 

one dollar will 
Mar 2, 1981 356 .00 

buy. 


(Israel's unit has been standardized to the early 1980 unit.) 


Clearly, both currencies have been inflating, at least as 
compared to the dollar. But the dollar itself has been 
inflating, at a rate over 10% in the same time period. 
There should be a way to calibrate a currency on a more 
solid basis than its value in dollars (perhaps against 

a weighted total of all the other currencies), and therein 
might lie a computing problem. 


5. Back in 1962 there was the PSAG project (the 
initials are those of Parkin, John Selfridge, George 
Armerding, and me). (Note: there are 24 possible 
permutations of the four letters, and we voted unanimously 
not to use either GASP or GAPS.) With advice and 
counsel from D. H. Lehmer, we had modestly set out to 
obsolete most of the number theory tables in the world. 
Specifically, we were going to calculate eight different 
funetions for each of the first 6,000,000 primes (for 
example, the least exponent, the least positive primitive 
root, and so on). A test run was made at UCLA, which 
indicated that our logic was correct and that Armerding's 
program was debugged and tested. Selfridge had devised 
a scheme whereby the complete factorization of (p - 1) 
for each prime p could be coded into 5 characters, which 
would be invaluable when it came time to print out the 
48,000,000 numbers we intended to produce. Parkin 
spent some happy hours with a punched card sorter, working 
out what we called Lehmer's Magic Numbers; they were the 
necessary conditions for small integers to be primitive 
roots, and they would save enormous amounts of computing 
time. A writeup on the PSAG project might have been 
interesting. To date, there has been no public outcry 
at the demise of the project. 


6. There has never been a thorough description of 
the Mendenhall-Warren-Hollerith method of progressive 
digiting; this was a 1929 scheme for developing sums of 
products using punched card gear of the time. The method 
was widely used from 1929 to around 1948, during which time 
it was continually improved. With the introduction of 
various calculating devices (particularly the IBM 604), 
digiting fell into disuse. But it returned with a roar 
in the 1950's, as one of the methods for converting 
Hollerith information on punched cards into binary, for 
card entry to computers. In the 1960's, on computers 
using card input, more digiting was performed each hour 
than had been done in the whole world from 1929 to 1948. 


This is typical of many things in computing; they 
i away for some time, only to reappear stronger than ever. 
It may be safe to say that we will never go back to cathode 
ray tube storage.) The digiting process will probably 
be back again. In any event, it is a fascinating foot- 
note to computing history, and it should be reported. 


7. There are books on compiler writing, and even 
college courses in the study of compilers, but where is 
there any solid literature on interpreters? If we have 
some 500,000 personal machines today, increasing currently 
at around 200,000 per year, each of which executes BASIC 
or Pascal interpretively, then it would seem that 
interpreters would rate considerable attention. If 
everyone meekly settles for the sloppy programming 
represented by the BASIC interpreters now running in the 
popular personal machines (the world was using properly 
written interpretive systems before the kids who wrote 
those interpreters were born), then we are going backwards. 
The people who know how such things should be done have 
been singularly quiet. 


8. The introduction of computers into cryptography, 
plus the recent attention given to public key systems, has 
directed all efforts toward substitution cipher systems. 
Little or no attention is paid to transposition ciphers, 
and they can be made quite secure. Moreover, they, too, 
lend themselves to computerization, and this matter should 
be looked into. 


9. Since we fear no man nor beast, there was also 
going to be a critical article on the "Computing" merit 
badge of the Boy Scouts. 


10. And, of course, there is a healthy backlog of 
new problems, solutions to old problems, and interesting 
correspondence. Perhaps some worthwhile articles may 
emerge from all this for the new POPULAR COMPUTING. 


By and large, it's been fun. Beginning about two 
years ago, it began to resemble work, which both the 
publisher and I find cuts into the fun. It is time once 
more to exit gracefully, while the audience is still 
applauding. 


